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Abstract 

Attractors of cooperative dynamical systems are particularly simple; for example, 
a nontrivial periodic orbit cannot be an attractor. This paper provides characterizations 
of attractors for the wider class of coherent systems, defined by the property that no 
directed feedback loops are negative. Several new results for cooperative systems are 
obtained in the process. 
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Introduction 



We consider differential equations 

cix 

-7- = F(x), x€X, t>0, (1) 
at 

where X c R" is convex, its interior is dense in X, and the vector field F : X ^MT extends 
to a vector field on an open set. The maximally defined solutions 1 1-> 0,(a), t >0,a e X 
generate the local semiflow O := {OJ^eR^. We refer to F (or {F,X,W), or (F,X,R",(^)) as 
a system. Dynamical notions are applied interchangeably to F and O. 

OF 

Many biological situations are modeled by cooperative systems: -^^ > Oif j i. The 
biological interpretation is that an increase of species / tends to increase the population 
growth rate of every other species j. In this case O is monotone, meaning it preserves the 
vector ordering. This causes the crude dynamics of a cooperative system to be compara- 
tively simple; for example, there are no attracting cycles and every orbit is nowhere dense 
(Hadeler & Glas [13], Hirsch [16]). 

Here we show that some of the dynamical advantages of cooperative systems extend 
to systems having a significantly weaker property: F is coherent (another name is positive 
feedback system) if whenever io,...,iy, v e { 1, . . . , n} are such that 

dFi, . 

iy = io, ik-i ^ ik and — ^ (1 < A: < v) 



then. 



dFi. . 

— — (x) does not change sign (1 <k<v) 
dxi^ 



and 

dF;. dFi , 

-^(x) ■ ■ ■ -^(x) > 0, (V;c G X). (2) 

Our chief combinatorial result. Theorem 10, shows that by permuting the variables x, 
and changing the signs of some of them, any coherent system can be transformed into a 
dynamically equivalent system (F,X, W, O) with the following properties: 

• F is not merely coherent, it has the stronger property of being quasicooperative: for 
any (/i, . . . , i„^ as above, each factor in the left hand side of (2) is > 

• if F is not cooperative, there exists a cooperative system (F^ , X^ , , O^), 1 <ni <n, 
such that the the natural projection 

n: r^R"', (xi,...,x„)h^(xi,...,x„,) 

maps X onto X^ and semiconjugates F to F^ and O to O^: 

n o F{x) = F^ o n(x), n o Of(x) = ^] o U(x) if 0,(x) is defined 
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Mild geometrical conditions on X guarantee that for each equilibrium p of F\ the restric- 
tion of F to := X n n~^(j!?) is equivalent to a quasicooperative system (Fp,Xp, R"""')- 
This is the basis for inductive proofs of our main results. 

We turn to our main topic, attractors. An attractor for F is a nonempty invariant con- 
tinuum A c X that uniformly attracts all points in some neighborhood of A. If the attraction 
is not necessarily uniform we talk instead of an attracting set. Three types of attractors A 
have received special attention: 

Point attractors: A is single point, necessarily an equilibrium. 

Periodic attractors: A is a cycle, i.e., a periodic orbit that is not an equilibrium. 

Strange attractors, often called "chaotic". This somewhat vague term signifies that A is 
neither an equilibrium nor a cycle, and usually that A is topologically transitive and exhibits 
"sensitive dependence on initial conditions". Some authors also require that periodic orbits 
be dense in A. 

This paper is motivated by the question: What kind of nonequilibrium attractors A can 
exist in coherent systems? Theorem 1 shows that A cannot be topologically transitive; 
Theorems 2 and 3 give further dynamical information. Other results apply to more general 
monotone local semiflows. 

Statement of results 

A set is finitely transitive for a system (or a local semiflow) if it is the union of the omega 
limit sets of finitely many of its points. 

Theorem 1 A finitely transitive attracting set A for a system (F, X, R") reduces to an equi- 
librium in the following cases: 

(i) F is coherent, and X is open in R" or relatively open in a coordinate half-space 

(ii) F is quasicooperative, and every point of A is strongly accessible in Xfrom above, or 

every point of A is strongly accessible in Xfrom below 

(iii) F is cooperative, and each point of A is strongly accessible in Xfrom above or below 

A stronger conclusion. Theorem 16, holds for cooperative systems. 

The following result requires no additional geometrical conditions on X: 

Theorem 2 If(F,X, W), n>lis a coherent system, every orbit is nowhere dense. 

Conjecture In a coherent system with n> \, every orbit closure has measure zero. Even 
for cooperative systems this is known only for n = l. 

An attractor is global if it attracts all points of X. An equilibrium is globally asymptot- 
ically stable if it is the global attractor. The following theorem needs X to be open: 
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Theorem 3 Let (F, X, R") be a coherent system with X open in W. Assume there exists 
a global attractor A. Then there exists an equilibrium, and if it is unique it is globally 
asymptotically stable. 

Proposition 14 extends a basic result previously known only for strongly order-preserving 
local semiflows. The development of the concept "attractor" is discussed in the Appendix. 

Motivations 

A coherent system is one whose interaction graph (defined below) has no directed negative 
loops. A more restrictive condition, for graphs that are not necessarily strongly connected, 
is the requirement that the graph has no wndirected negative loops: in that case, one may 
always perform an elementary change of variables (defined below) that transforms such 
a system into a cooperative one. In a classical and often-quoted 1981 paper, R. Thomas 
conjectured that coherent systems do not have any periodic attractors: "the presence of at 
least one negative loop in the logical structure appears as a necessary (but not sufficient) 
condition for a permanent periodic behavior" [48]. It has often been claimed (see e.g. [30]) 
that Thomas' conjecture was settled in [41, 12]. However, these references only dealt with 
the more restricted monotone case. Theorem 1 in this paper settles the question. We refer 
the reader to [42] for further comments on the relevance of these concepts to molecular 
systems biology, and to [43] for numerical simulations which suggest that systems that 
are "close" to having the coherence property might have, in some statistical sense, simpler 
attractors. 

Structure of proofs 

The proofs of Theorems 1, 2 and 3 have a common pattern which we now discuss. Let T 
stand for one of these theorems. It is proved first for a cooperative system, which includes 
the case n = I. The proof proceeds by induction on n. A coherent system which is not 
cooperative is transformed, by permuting and changing signs of variables, to a system 
(F, X, R") having the following properties: 

• F is quasicooperative 

• there is a system (F\X\ R"') with ni < n, such that the natural projection 11: R" — > 
R"! satisfies 

U(X) =X\ no F(x) = Fo n(x), (x e X) (3) 

• FMs cooperative 

It follows that n semiconjugates the local semiflow O of F to the local semiflow of F^: 

n o 0,(x) = o n(x) if OXx) is defined (4) 
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We summarize this by saying that U: (F,X,R") (F\X\R"') (or U: F ^ F^) is a 
cascade. We also allow the trivial cascade, for which F = F^. 

For each equilibrium p of F^ the affine subspace Ep := Il~^(p) is a coset of the kernel 
of n. The canonical chart 

Tp'. Epx'R""\ (xi, . . . ,x„) ^ (x„^+\, . . . ,x„) (5) 

is an affine automorphism. 

The vector field F, being tangent to Ep along Xp, restricts to a vector field Fp inXp := 
Xp :=Xr\Ep, and O restricts to a local semiflow in X^. The hypothesis of T will ensure 
that the relative interior of Xp in Ep is dense in Xp. The canonical chart converts the fibre 
system (Fp, Xp, Ep) into a system (Fp, Xp, R."""')- 

We identify each fibre system (Fp, Xp, Ep) with (F, Xp, R"""') by means of the canonical 
chart. Thus Fp has an interaction graph Y(Fp) := T(Fp). We ascribe to Fp the property of 
being cooperative, quasicooperative or coherent whenever that property holds for Fp. 

Theorem 7~ holds for the cooperative system F\ and it holds for all fibre systems by 
the inductive assumption. The induction is completed by showing that this implies 7~ also 
holds for (F,Z,R"). 

There is a delicate point regarding the domains of these systems. The proofs for coop- 
erative systems use special properties of X, such as every point being strongly accessible 
from above. These properties are postulated in the hypotheses of the main theorems. To 
make the induction work, the same properties must be verified for the systems obtained by 
elementary coordinate changes, and also for fibre systems. This means that the class of 
domains X referred to in the theorems must be preserved by permuting and changing signs 
of variables, and by intersecting X with the affine subspaces Ep. For this reason X is usually 
required to be an open set in R" or a relatively open subset of a coordinate halfspace. 

Local semiflows 

A local semiflow O in a metrizable space Z is a collection O = {^t}t€K+ of continuous maps 
Of : Dt ^ Rt between nonempty subsets of Z, with Dt open. The notation 0,x indicates 
jc G Dt, absent contraindications. O is required to have the following properties: 

• The set := {(t, x) G R+ x Z : x G D,} is an open neighborhood of {0} x Z in R+ x Z, 
and the map Q.^ Z, (t, x) i-> 0,x is continuous. 

• X G (<D,)-iZ), => o 0,(x) = 0,«(x) 

• Oo is the identity map of Z. 

We also say that (0,Z) is a local semiflow. When O is obtained by solving Equation (1) 
each map is a homeomorphism, but this is not assumed for general local semiflows. 
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The orbit and omega limit set of x are respectively 

y(x) := {Q>tx: x e D,], oj(x) := p| r(0,^) 

t>Q 

p is an equilibrium if O,/? = /? for all t. The set of equilibria is denoted by £(0), and by 
£>{F) when O is generated by the vector field F. 

Attractors and attracting sets 

We call A positively invariant for if 0,((2) is defined and belongs to A for all ? > 0, a e A, 
and invariant if in addition A is nonempty and Of(A) = A for all t > Q. We say that A 
attracts x if 7(jc) is compact and c A. The set of such points y is the basin of A. 

A is topologically transitive if it is the omega limit set of one of its points, and finitely 
transitive if it is the union of the omega limit sets of finitely many of its points. 

We call A attracting if it is invariant, connected and compact, and its basin is a neigh- 
borhood of A. If in addition A has arbitrarily small positively invariant neighborhoods, A is 
an attractor} 

Ordered spaces 

By an ordered space we mean a topological space Z together an order relation K dZxZ 
that is topologically closed. If x, y e Z we write: 

x>y and y < x if {x,y) e H, x>y and y<x if x>y,xi^y (6) 

The vector order in any subspace of R" is defined by 

u>v <==^ u-veKl 

where R" denotes the the positive orthant [0, oo)" c R". 

A subset of an ordered space is unordered if none of its points are related by >. 

Every subspace X <z Z inherits an order relation from Z. If M c Z then x > M means 
x> yfox all y & M, and similarly for the other relations in (6). For jc, j e X we write 

x>xyif x> N, y e Intx(AO, 
x<xyif X <N, y e Intx(AO 

for some open N cX. Note the notational anomaly that x >x y and y<xxme not equivalent 
statements for general ordered spaces. They are equivalent, however, if Z c R" is open and 
has the vector ordering. For example, in Z = R^ we have (0, 0) <x (0, 1) but (0, 1) />x(0, 0). 

'There are many definitions of "attractor" in current use, not mutually consistent. The one adopted here 
is equivalent to that of Conley [10], and (for compact invariant sets) those of Hale [14] and Sell & You [38]. 
It is analogous to the definitions for discrete-time systems in Smale [39] and Akin [1]. 
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Let X be a subset of an ordered space Z. We call q € X strongly accessible in Xfrom 
above (respectively, from below) if every neighborhood of ^ in Z contains a point x>x q 
(respectively, x <x q)-^ 

All our results are valid when X is an open set in R", and some are valid for special kinds 
of nonopen sets, especially open subsets of a coordinate halfspace of R", which means a 
set 

{x G R" : axi > q} 

for some choice ofl&{l,...,n},ae {+1}, (ci, . . . , c„) G R". We rely on the following fact, 
whose proof is left to the reader: 

Lemma 4 Assume X c R" has the vector ordering. If X is an open subset of R", or a 
relatively open subset of a coordinate halfspace, every point ofX is strongly accessible 
from above and below in X. | 

Note also that if X is an open subset of R" , all points of X are strongly accessible in X from 
above. 



Cascades 

Let (F,X,R") and (i^\X\R"') be systems with I < ni < n and assume 11: F -» FMs a 
cascade (see (3)). This implies 



— = if i<ni< ], 
dxj 



(i,7 e {1, 



(7) 



and the Jacobian matrices of F have lower triangular block decompositions of the form 



F'(x) = 



Mn(x) O 
M2i(x) M22(X) 



(8) 



where Mii(x) = (F^YiUx) e R"'^"', and O stands for a matrix of zeroes. The following 
diagrams commute for each r > 0: 



Dm 

n 

Z)(Oi) 



n 



For p G S(F^) let Tp-. Ep x R"'"' be the canonical chart. Set Xp = Tp{Xp) and define 
Fp-. Xp ^ R"~"i to be the unique vector field transformed by (Tp)~^ to Fp, that is. 



^Slightly stronger properties with the same names are used in Hirsch & Smith [20]. 
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F{TpX) = TpF(x), (x G Xp) (9) 

The local semiflows of Fp and Fp are conjugate under Tp. For (Fp,Xp, R"~"i ) to be a system 

it is necessary and sufficient that the relative interior of Xp in Ep be dense in Ep. When this 
holds we call (Fp, Xp, Ep) the fibre system over p and identify it with (Fp, Xp, R" ) by Tp, 

The interaction graph T(Fp) := T(Fp) is determined by the signs of the entries in the 
block M22(x) in (8). The next lemma gives convenient conditions ensuring this. 

Consider the following conditions: 

• Ci(X,R"):XisopeninR" 

• C2(X, R"): X is open in a coordinate half space of R" 

• C3(X,R"):ZisopeninR:^ 

• C4iX, R"): Z is a rectangle 

Lemma 5 Assume 11: F ^ F^ a cascade as above and p e 8(F^). Suppose p e S(F^), 
and CdiX,R"), is satisfied for some d e {1,2,3,4}. Then (Fp,Xp,R"~"^) is a system, and 
Cd(Xp,R"-''') holds. 

Proof The verification that C^iX, R") impUes C^iXp, R"-"0> and also that relative interior 
of Xp in Ep is dense in Xp, is straightforward. I 



Graphs 

By a directed graph Y := (Vr,Er) we mean a nonempty finite set V := Vr (the set of 
vertices) together with a binary relation E := Er <z V xV (the set of directed edges, usually 
referred to simply as "edges"). We always assume E is totally nonreflexive i.e., (i, i) i E. 

An isomorphism between a pair of directed graphs is a bijection / between their vertex 
sets such that / x / restricts to a bijection /, between their edge sets. 

Our chief tool for analyzing the crude dynamics of systems (F, X, R") is the interaction 
graph r := Y{F). This is the labeled directed graph with vertex set is V = V(r) := 
{\,...,n], whose set of (directed) edges is 

dFi 

E = E(r) := {(j, i)€VxV: 7 / and — ^ is not identically in X} 

dxj 

Edge (7, is assigned the label h(j, i) e {+1,-1,6} according to the rule: 



Kj, i) = 



1 if Tiix) > for all x £ X, 



dx-i 



-1 if ^(x) < for all x£X, (10) 



dxJ ' 

6 otherwise 
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and is respectively called positive, negative or ambiguous. A loop is positive if each of its 
edges is labeled +1 or -1 and the product of these labels is +1. 
We define three types of graphs in increasing order of generality: 

r is positive if every edge is positive, 

r is quasipositive if every loop has only positive edges, 

r has the positive loop property if every loop is positive. 

Paraphrasing some of the earlier definitions, we define corresponding types of systems F 
in terms of T{F): 

F is cooperative if r(F) is positive 
F is quasicooperative if T{F) is quasipositive, 
F is coherent if Y has the positive loop property 
Evidently cooperative => quasicooperative => coherent. 

The term "graph" is shorthand for "finite directed graph having edges labeled in {1, -1,0}." 
Graphs are denoted by Greek capitals F, A, perhaps with indices. The sets of vertices and 
edges of T are denoted by V(Y) and ^(r), respectively, and the labeling function is denoted 
by hi : V(Y) — > {1,-1,6}. Two graphs Y, A are isomorphic if there there is an isomorphism 
/ : V{Y) —> V(A) between the underlying directed graphs such that /ia o = hr. 

A is a subgraph of Y provided 

y(A) c V(Y), E(A) c E(Y), hj, = hr\E(A), 

We abuse notation and denote this by A c F, saying that A contained in F. 

If A, A' are subgraphs their graph union is the subgraph with vertex set V(A) U V(A') 
and edge set E(A) U E(A'). 

A path of length ^ e N+ is a sequence (uq,..., Uk) of vertices such that (v^-i, v^) is an 
edge for j = l,...,k. The concatenation of an ordered pair (A,p) of paths, 

A = (Uq, . . . , Uic), p = (Uic, . . . , Uk+i), 

is the path 

A- p := (uo, . . .,Uk, Uk+i, . . . , Uk+i) 

obtained by transversing first A and then p. 

A loop of length // e N+ is a sequence of p>2 edges having the form 

Oo, ii), 0"i, ii), (i/x-i, i/x), i/x = io 

As our graphs are totally nonreflexive, there are no self-loops: ij ^ i, j - I,..., p. 

A loop is positive (respectively, negative) if each of its edges is labeled 1 or -1 and the 
product of these labels is +1 (respectively, -1). All other loops are ambiguous. 

In the next three definitions the labeling plays no role. A graph is called: 

connected if for each pair of distinct vertices j, k there is a sequence of vertices 
j = io,. . ., im = k, m & N+ such that (j/-i, i/) or (j/, is an edge of A, (/ = 1, . . . , m) 
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strongly connected if for any ordered pair (a, b) of distinct vertices there is a path 
in A from a to b, 

primary if every edge belongs to a loop. 
These definitions imply: 

• A graph with no edges is primary, but a graph with only one edge is not primary. 

• The graph union of primary subgraphs is primary. 

• A strongly connected subgraph is primary, and a primary connected subgraph having 
more than one vertex is strongly connected. If Y is quasipositive, every primary 
subgraph is positive. 

A subgraph A c F is called: 

• full provided it contains all edges in F joining vertices of A, 

• initial if no edge of F is directed from a vertex outside A to a vertex of A, 

• terminal if if no directed edge of F joins a vertex of A to a vertex not in A, 

• fundamental if is connected, primary and initial, and no other subgraph containing 
A has these properties. 

Lemma 6 The following hold for all subgraphs: 

(a) fundamental subgraphs are full 

(b) if fundamental subgraphs Ai, A2 share a vertex, they coincide 

(c) every connected, primary, initial subgraph is contained in a unique fundamental sub- 

graph 

Proof (a) and (b) follow directly from definitions, (c) is proved by showing that the graph 
union of a maximal nested family of connected, primary, initial subgraphs is fundamen- 
tal. I 

Graphs and systems 

Let (F, X, R") be a system. 

Proposition 7 If II: F ^ is a cascade having a fibre system Fp, then: 

(a) T(F^) is a full subgraph ofY{F). 

(b) T{Fp) is isomorphic to a subgraph ofT{F). 
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(c) when F is cooperative, quasicooperative or coherent, F and Fp have the same prop- 
erty. 

Proof (a) and (b), which imply (c), are proved by inspecting the block decomposition (8) 
of the matrix of functions F'(x). I 

Propositions Let c T(F) be an initial full subgraph such that V(r{) = {!,..., ni}. 
Then: 

(i) there is a cascade U : (F, X, W) (F\X\R"') such that r(F^) = T\ 

(ii) When F is quasicooperative, F^ and all fibre systems are quasicooperative, and ifYi 

is primary then F^ is cooperative 

Proof Initiality and fullness of means that (7) holds. Therefore (3) defines a cascade 
satisfying (i). The first assertion in (ii) follows from Proposition 7(c). The second assertion 
holds because Fi is quasipositive, and if it is primary all its edges are in loops and hence 
are positive. I 



Spin assignments 

A spin assignment for a graph F is any function cr: y(F) — > {+1}. It is consistent if 
h(u,v) = cr(u)cr(v) for every edge (u,v) belonging to a loop. (This terminology is not the 
same as in [42], where it was required that every edge be consistent. With that stronger 
requirement, the theorem given below would become a characterization of monotonicity 
with respect to an orthant order, a more restrictive property than coherence.) 

Theorem 9 F has the positive loop property if and only if it has a consistent spin assign- 
ment. 

Proof Assume F has the positive loop property. Let F' be obtained from F by keeping the 
same vertices but deleting the edges not contained in loops. Clearly F' has the positive loop 
property, and if cr is a consistent spin assignment on F' it is also consistent on F. Therefore 

we can assume every edge e belongs to a loop and is thus positive. 

Claim: If A^,A^ are paths from aiob then h{A^) = h{A^) e {±1}. To see this, choose 
a path fi from b to a, which can be done because each edge belongs to a loop. Since every 
loop is positive by hypothesis, for j = 1,2 we have 

1 = hiAJ ■ lu) = hiAj)h(ju) 

Therefore h(A^) = h(ju) = h(A^). 

Now fix a vertex p of F and for each vertex v choose a path Ay from p to v. Define 
<t(p) = 1 and cr(v) = h(Ay), which by the claim is independent of the choice of Ay. For any 
edge e = (u, v) we can fix Au and define Ay := Au- e. Then have: 

o-(m) = h(Au), cr(v) = h(Au ■ e) = h(Au)h(e), 

which implies h(e) = o-(u)cr(v). The converse implication is left to the reader. | 
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Remark The foregoing proof can be expressed homologically. Let A denote the 1- 

dimensional cell complex corresponding to a prime subgraph A c F having the vertices 
of A for 0-cells and the directed edges of A for 1 -cells. In the cellular chain groups of 
A with coefficients in Z2 (identified with the multiplicative group {±1}), a labeling h is a 
1-cochain, spin assignments are 0-cocycles, and a spin assignment cr is consistent for h if 
its coboundary is 60- = h. As the evaluation of cochains on chains induces a dual pairing 
H^(A; Z2) X Hi(A; Z2) — > Z2, the positive loop property makes the cohomology class of h 
trivial. Thus h = 5a, proving that cr is consistent. 

A change of variables x y h called elementary if there is a permutation / 1-^ /' of 
{\,...,n} and an n-tuple p e {+1}" such that j,- = p;jc,v. 

Theorem 10 Tjf a system is coherent, there is an elementary change of variables trans- 
forming it to a quasicooperative system admitting a cascade over a cooperative system for 
which all fibre systems are quasicooperative. 

Proof Assume (F, X, R") is a coherent system, which by Theorem 9 has a consistent spin 
assignment cr. The elementary change of variables L: R" ^ R", 

y = Lx, yi := cr(i)xi 

transforms (F, X, R") into a system 

(G,L(X),R"), LoG*** 

such that r(G) and T{F) have the same undirected edges. For every directed edge (7, /) of 
r(G): 

^r(G)(j, = sign(-^) = o-^CT; sign(-^) crjO-ihr(F)U, 
If (j, i) belongs to a loop then hr{F)(j, i) = crjcri by the consistency condition. Therefore 

sign(^)((r,.cr,)2 = (+1)2 = 1, 

showing that G is quasicooperative. After reindexing variables we assume there is a funda- 
mental subgraph F^ c F(F) with vertex set {!,..., ni), I < ni < n. Now apply Proposition 
8. I 

Monotone dynamics 

A local semiflow O is monotone if x > y => > O^y. Throughout this section we 
assume: 

• O := {Of}f>o is a monotone local semiflow in an ordered space X 
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To simplify notation we may write x(t) := 0,x whenever O^x is defined. It is well known 
for the data in Equation (1) that if F is cooperative and X is convex, the corresponding 
local semiflow O is monotone. This is a corollary of the MuUer-Kamke theorem [28, 21] 
on differential inequalities (Hirsch [16]). 

Proposition 11 The following are true for all x G X: 

(a) No points ofcoix) are related by >x or <x 

(b) cij(x) is a singleton in the following cases: 

(i) 7(x) is compact and there exist > 0, e > such that 

t^ < t < t^ + s => 0,;c < xor <D;;c > x 



(ii) y(x) is compact and there exist t > such that 

(t>tX >x X or ^fX <x X 

Proof (a) and (b)(i) are sharpenings of Hirsch & Smith [20, Theorems 1.8, 1.4], respec- 
tively. Assertion (b)(ii) follows from (b)(i). I 

Proposition 12 Assume A c X is attracting. 

(a) If each point of A is strongly accessible in Xfrom either above or below, then A contains 

an equilibrium. 

(b) If each point of A is strongly accessible in Xfrom both above and below and A OS = p 

then A = p. 

Proof This is a slight generalization of Hirsch [17, Theorems ni.3.1 and in.3.3], and the 
same proofs work here. I 

Proposition 13 Assume A c (i>(x). Let q ^ Abe a minimal (respectively, maximal) point 
of A having a neighborhood N cX such that there is a point y < N ( respectively, y > N) is 
attracted to A. Then q = inf A (respectively, q = sup A). 

Proof To fix ideas we assume ^ is a minimal point of A and y < N. Notation is simplified 
by setting O^w = w(t) whenevever w e X,t > 0. 

Some point on y(x) lies in N its omega limit set contains A. Replacing x by such a point 
we assume x e N. Therefore y < x and 

y(t)<x(t), (t>0) (11) 
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There is a sequence ^ oo such that x(?„) e and 

x{t„) (12) 

Because (o(y) meets A we can choose this sequence so that also 

y{tn) ^ fl G A (13) 

It follows from (11), (12), (13) and closedness of the order relation that a < ^, so minimality 
of q implies a = q. Thus 

y{tn)^q&& (14) 

Choose hq so that y{tn^) G. N. If / c R+ is a sufficiently small open interval about t„g 
then s e I => y(s) e N, hence y(s) > y. The dual of Proposition ll(b)(i) now shows 
that a)(y) is an equilibrium, hence oj(y) = {q}. It follows from (11) that oj(x) > q, hence 
A>q. I 

In the rest of this section we assume: 



• X c R" with the vector ordering. 

Proposition 14 Assume x e a>(A), a>(x) = A. If inf A = p or sup A = p then A = p. 

This result also holds when X is ordered by a solid polyhedral cone, but it is has not been 
proved for more general ordered spaces. For strongly order-preserving local semiflows 
a stronger conclusion holds: Every omega limit set is unordered (Hirsch & Smith [20, 
Corollary 1.9]). 

Proof For any Z c {1, . . . , n} the corresponding/ace of R" is 

J := J(S) = {z e r:^ : Zi>0 => / e E} 
When E the corresponding open face is 

J" := J°(i:) = {z e r: : z,- > « i e 1} 
It can be seen that J" = J and J" is relatively open in its linear span. Moreover 

(Vz e J") (35 > 0) z>JnA^5(0) (15) 

Fix x,peX such that 'mia){x) = p or sup cl>(x) = p; we have to prove oj{x) = p. To 
fix ideas we assume p = = mftoix). Claim: 0,(.)c) is defined for all t > 0. It is well 
known that this is the case if the orbit closure of x is compact. If it is not compact, the orbit 
intersects the boundary of some open ball centered at in an infinite set. Consequently 
a>(x) contains a point p, which implies the claim. 

For any / c [0, oo) set <!>(/, x) := {0,jc: t e I}. By the Baire category theorem there is a 
dense open subset S c [0, oo) such that for each component I of S there is a unique open 
face J^^ D <!>(I,x). 
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There is a sequence {4} of these components and points tk g h such that as ^ ^ oo we 
have 

tk oo, x(tk) > 0, x(tk) 

After passing to a subsequence we can assume there is an open face K" such that J^^ = K" 
for all k. Choose such a K" having the largest possible dimension. Then x{t) e K" for 
sufficiently large t. For if x(to) e K'' and e > is such that x(t) i K" for t e (?o, t^-^ s\, 
then x{t') for some t' G (fo, h + ^] belongs to an open face of larger dimension, and this can 
only happen finitely many times. Set dim K" = m G {1, . . . , n} and relabel variables so that 
K'' = K''({l,...,m}). 

By (15) there exists > such that 

? > 4 => x{Q > x(t) > 

By Proposition 1 l(b)(i) the trajectory of x(U) converges, necessarily to 0. Therefore <jj(x) = 
a>(x(Q) = 0. I 

Corollary 15 Assume eX cRl. If e a>(x), then = a>(x). I 

Proof Follows from Theorem 14 because = inf a)(x). I 



Remark We digress to interpret this result biologically. Let Xi > stand for the "size" of 
species / (population, biomass, density, . . . ) and call ^Li "total size". Assume that 
from each initial state jc(0) g R" the species develop along a curve x(t) = (xi(t),. . ., x„(0) g 
R+,^ ^ governed by a cooperative system (suggesting symbiosis or commensalism) in 
R!^. Then: 

• If the total population does not die out, the total size is bounded above 0. 

This follows from the contrapositive of the Corollary. 

The next result will be used to start the inductive proof of Theorem 1. It applies only 
to cooperative systems, but the assumptions on O, X and A are weaker than in Theorem 1. 
Recall that every nonempty compact set in an ordered space contains a maximum point and 
a minimum point (Ward [51]), 

Theorem 16 Assune X c R" has the vector ordering and O is a monotone local semiflow 
in X. Let A c X be attracting and finitely transitive for O. If every point of A is strongly 
accessible in Xfrom above or below, then A € S. 

More precisely: If q € A is maximal and strongly accessible in X from above then 
A = q. Likewise ifq€A is minimal and strongly accessible in Xfrom below. 
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Proof It suffices to assume q £ A is maximal and strongly accessible in X from above. 

Under the current assumptions there exist x £ A,y e X and neighborhood U c X of q such 
that that q G co(x), y > q and y is attracted to A. Evidently q is maximal in a)(x), hence 
q = sup cl)(x) by Proposition 13, and therefore q = a>(x) by Theorem 14. 

Suppose z e X and oj(z) (MJ <Z. There exists / e N+ with zih) e U, hence y{t) >xz{t + 
0, {t > 0), and monotonicity proves 

oj{z) nU ^ coiz) < q (16) 

Now we prove for all v e X: 

q e a>(v) => q = co(v) (17) 

For there exists z G 7(v)n ?7 and Equation (16) implies q = sup <x>(z), hence q = a>(z) = cl>(v) 
by Theorem 14. 

Let {a^:} be any sequence in U n A converging to q. By hypothesis there is a finite set 
S a X such that each aj^ is an omega limit point of some member of 5 . By finiteness of S 
there is a subsequence {bk} of {ak) and v G. S such that {^^t} c a>(v). Evidently q G a»(v), 
whence q = (X»(v) by (17). This can only happen ifbk = q for all A: G N+. It follows that q is 
isolated in the connected set A, entailing A = q. | 



Proofs of the main theorems 

Proof of Theorem 1 Let the system (F, X, R") be as in Theorem 1, with a finitely transi- 
tive attracting set A <zX. 

Step ( i) Consider first the case that F is cooperative. Then O is monotone because X is 
convex, and each of the assumptions (i), (ii) implies each point of A is strongly accessible 
in X from above or below. The conclusion for this case follows from Proposition 16. 

Step (ii) We proceed by induction on n, the case n = I following from the cooperative 
case. Assume inductively that n > I and that the conclusion holds for smaller values of 
n. By Step (i) we can assume F is not cooperative, whence by Theorem 10 there is a 
cooperative system (F\Z\R"') and a cascade II: F ^ F^ with 1 < ui < n, whose fibre 
systems are quasicooperative. Lemma 5 shows that (Fp,Xp, Bp) is a fibre system for each 

The set n(A) c X^ is finitely transitive for the cooperative system F^ , hence n(A) = 
p G £(F^) by Step (i). Thus A lies in the invariant set Xp =Xf\ Q'^iA), and A is attracting 
and finitely transitive for := Q>\Xp. The inductive hypothesis applied to (Fp,Xp,Ep) 
shows that A is an equilibrium, completing the induction. | 

Proof of Theorem 2 Consider first the case that F is cooperative. Assume per contra 
that the orbit closure of x e X contains a nonempty open subset U c X. As some open 
subset of W is dense in X we can assume U is open in R'\ The orbit y(x), being a smooth 
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curve, is nowhere dense in U because n > 2. Therefore U c a>(x), hence a)(x) contains 
points a, b such that a >x b. But this contradicts Proposition 1 1(a). 

Now assume F is not cooperative. By Theorem 10 there is a cascade 11: F ^ with 
F^ cooperative. liW <zX\s open and y is an orbit of F, then n(W) is open in and 11(7) 
is an orbit of The cooperative case shows that 11(7) Pi Il{W) is not dense in Il{W) and 
therefore y n W is not dense in W. I 



Proof of Theorem 3 If F is cooperative, as when n = 1, the conclusion follows from 
Proposition 12. We proceed by induction on n, assuming that n > I and the theorem holds 
for smaller values. 

We can assume F is not cooperative. By Theorem 10 there is a cascade 11: F ^ F^ 
cooperative system (F^,X^,R"^) with F^ cooperative and 1 < «i <n, such that if p e £(F^) 
then (Fp,Xp,Ep) is a quasicooperative system. Applying the inductive hypothesis twice, 
we conclude that there exists p e G{F^) and q e &{Fp) c &{F). 

Assume &{F) = q and set n(^) = p' & £>{F^)- Then p' = p. For we showed above 
that every fibre system contains an equilibrium of F, which must be q. Thus ir\p') and 
Il~^(p) are not disjoint, hence they coincide and 11 maps both of them to p. 

By the inductive hypothesis p is the global attractor for F\ therefore Xp attracts all 
points of X by Equation (4). This implies A is the global attractor for 0|Xp, and the induc- 
tive hypothesis applied to (Fp, Xp, Bp) shows that A = q. I 



Appendix: Notes on the development of the concept "at- 
tractor" 

In spite of the fact that everyone who is interested in dynamics has a more or less vague 
intuition of what an attractor of a map f:M—>M should be, there is no generally 
accepted mathematical definition for this concept even if M is a smooth manifold and 
f is also smooth. — H. Bothe [4] 

The first mathematical use of the word "attractor" may be in Coddington & Levinson's 
1955 book [9], where it refers to an asymptotically stable equilibrium. The term was subse- 
quently extended to include an attracting cycles. Today there are many definitions, usually 
meaning an invariant set (of some kind) that is approached uniformly (in some sense) by 
the forward orbits of all (or most) points in some neighborhood of the set. 

Attractors do not occur explicitly in the work of Poincare or Birkhoff". These authors 
were primarily interested in Hamiltonian systems, which have no attractors because they 
preserve volume. 

An early proof of existence of a unique attracting periodic orbit for a general class of 
systems is in the 1942 paper of N. Levinson and O. Smith [23]. ^ 

^Thanks to George Sell for this reference. 
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Early computer simulations revealed what appear to be attractors. As far back as 1952, 
Turing [49] published pictures of numerical simulations of a nonlinear dynamical model 
of cell development, exhibiting striking pattern formation. Simulations by Stein & Ulam 
[44, 45] and Lorenz [24] gave persuasive pictorial evidence of complicated structure in 
attractors, but attracted little attention when they were published. Hamming's review [15] 
of [45] was unenthusiastic: 

Many photographs of cathode ray tube displays are given, a fondness for citing large 
numbers of iterations and machine time used is revealed, and a crude classification of 
the hmited resuhs is offered, but there appears to be no firm new results of general 
mathematical interest. . . 

One can only wonder what will happen to mathematics if we allow the undigested 
outputs of computers to fill our literature. The present paper shows only slight traces 
of any digestion of the computer output. 

Much of the early theoretical work on attractors on global analysis was concerned with 
characterizing them in terms of Liapunov functions and topological dynamics (e.g., Ura 
[50], Auslander et al. [2], Mendelson [25], Bhatia [3]). Little was known of their internal 
dynamics beyond the existence of fixed points in global attractors for flows in Euclidean 
space (Bhatia & Szego [5]). 

In the 1960s a number of articles on attractors and related forms of stability were in- 
spired by Sell [37]. In his seminal 1967 work on global analysis, Smale gave detailed 
constructions and analyses of hyperbolic attractors and other invariant sets, which would 
later be called "chaotic" and "fractal", and proved them structurally stable. He called at- 
tention to the vast mixture of periodic, almost periodic, homoclinic and other phenomena 
found in structurally stable attractors, even in rather simply given systems. 

"Strange attractors" were proposed in 1971 as a model of turbulence by Ruelle and 
Takens [35, 36, 32], Newhouse et al. [29]). The physical significance of this route to chaos 
is still debated. 

In his controversial 1972 book on morphogenesis ([46, 47]) the late Rene Thom issued 
a bold manifesto proclaiming the fundamental scientific role of attractors: 

1. Every object, or every physical form, can be represented by an attractor C of a 
dynamical system in a space M of internal variables. 

2. Such an object possesses no stability, and for this reason cannot be perceived, unless 
the corresponding attractor is structurally stable. 

3. Every creation or destruction of forms, every morphogenesis, can be described by 
the disappearance of the attractors representing the initial forms and their replacement 
through capture by the attractors representing the final forms. This process, called 
catastrophe, can be described in a space of external variables. . . . 

In recent years much work has been devoted to analysis of attractors in specific classes 
of chaotic systems, such as those named after Duffing, Lorenz, Henon and Chua, and to at- 
tractors having particular topological properties, such as R. Williams' expanding attractors 
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(Williams [52], Plykin & Zhirov [31]). A novel measure-theoretic type of attractor due to 
Milnor [27] has stimulated several papers. 

Many authors have investigated attractors in infinite-dimensional systems, especially 
for partial differential equations, a prime desideratum being finite dimensional gUobal at- 
tractors. The large literature includes books by Constantin etal. [11], Hale [14], Ladyzhen- 
skaya [22], Ruelle [33], Sell & You [38], and others. 

Attractors, being objects defined by topological limiting processes, resist classification 
and even description. A general theory appears quite distant. 

References 

[1] E. Akin, "The general topology of dynamical systems." Providence: American Math- 
ematical Society (1993) 

[2] J. Auslander, N.P Bhatia & P. Seibert, Attractors in dynamical systems, Boletin Soc. 
Mat. Mex. (2) 9 (1964), 55-66 

[3] N. Bhatia, On asymptotic stability in dynamical systems. Math. Systems. Theory 1 
(1967), 113-128 

[4] H. Bothe, Review of Milnor [27], Math. Reviews. 

[5] N. Bhatia & G. Szego, Weak attractors in Math. Systems. Theory 1 (1967), 129- 
133 

[6] V. Chepyzhov & M. Vishik. "Attractors for equations of mathematical physics." 
American Mathematical Society Colloquium Publications, 49. Providence, RI: Amer- 
ican Mathematical Society (2002) 

[7] J. Cholewa & T. Dlotko, "Global attractors in abstract parabolic problems." London 
Mathematical Society Lecture Note Series, 278. Cambridge, England: Cambridge 
University Press (2000) 

[8] S. N. Chow & J. K. Hale, "Methods of bifurcation theory." New York: Springer- Verlag 
(1982) 

[9] E.A. Coddington & N. Levinson, Theory of ordinary diff'erential equations. New 
York: McGraw Hill (1955). 

[10] C.C. Conley, "Isolated Invariant Sets and the Morse Index." Published for the Con- 
ference Board of the Mathematical Sciences by the American Mathematical Society, 
Regional conference series in mathematics. No. 38, Providence (1978) 



19 



[11] p. Constantin, C. Foias, B. Nicolaenko, R. Temam, "Integral manifolds and inertial 
manifolds for dissipative partial differential equations." Applied Mathematical Sci- 
ences, 70. New York: Springer- Verlag, 1989. 

[12] J.L. Gouze, Positive and negative circuits in dynamical systems, J. Biol. Sys., 6 (1998) 
11-15. 

[13] K.P Hadeler & D. Glas, Quasimonotone systems and convergence to equilibrium in 
a population genetic model, J. Math. Anal. Appl. 95 (1983), 297-303 

[14] J.K. Hale, "Asymptotic behavior of dissipative systems." Math. Surveys and Mono- 
graphs 25. Providence, R.I.: Amer. Math. Soc. (1988) 

[15] R. Hamming, Review of Stein & Ulam [45]. Math. Reviews # 6666 29 (1965), 1248. 

[16] M.W. Hirsch, Systems of differential equations which are competitive or cooperative. 
I: limit sets. SIAM J. Math. Anal. 13 (1982), 167-179 

[17] M.W. Hirsch, The dynamical systems approach to differential equations. Bull. Amer. 
Math. Soc. 11 (1984), 1-64 

[18] M.W. Hirsch, Stability and convergence in strongly monotone dynamical systems, J. 
reine und angewandte Mathematik 383 (1988), 1-53 

[19] M.W. Hirsch, Chain transitive sets for smooth strongly monotone dynamical systems, 
Dyn. Contin. Discrete Impulsive Systems, 5 (1999), No. 1-4, 529-543 

[20] M.W. Hirsch & H.L. Smith, Monotone Dynamical Systems, "Handbook of Differen- 
tial Equations: Ordinary Differential Equations, Vol. 2." A. Canada, P. Drabek & A. 
Fonda editors, 239-258. Boston: Elsevier North Holland (2005) 

[21] E. Kamke, Zur Theorie der Systeme gewohnlicher Differentialgliechungen 11, Acta 
Math. 58 (1932), 57-85. 

[22] O. Ladyzhenskaya, "Attractors for semigroups and evolution equations," Cambridge, 
England: Cambridge University Press (1991) 

[23] N. Levinson & O. Smith, A general equation for relaxation oscillations, Duke Math. 
J. 9 (1942), 382^03 

[24] E. Lorenz, Deterministic nonperiodic flows, J. Atmos. Sci. 20 (1963), 130-141 

[25] P Mendelson, On unstable attractors, Boletin Soc. Mat. Mex. (2) 5 (1960), 270-276 

[26] J. Mierczynski, A remark on M. W. Hirsch's paper "Chain transitive sets for smooth 
strongly monotone dynamical systems," Dynam. Contin. Discrete Impuls. Systems 7 
(2000), 455-461. 



20 



[27] J. Milnor, On the concept of attractor. Comm. Math. Phys. 99 (1985), 177-195 

[28] M. Miiller, Uber das fundamenthaltheorem in der theorie der gewohnlichen differen- 
tialgleichungen. Math. Zeit. 26 (1926), 619-645. 

[29] S. Newhouse, D. Ruelle & F. Takens, Occurrence of strange Axiom A attractors near 
quasiperiodic flows on T"',m > 3, Comm. Math. Phys. 64 (1978-79), 35-40. 

[30] S. Pigolotti, S. Krishna, and M.H. Jensen, Oscillation patterns in negative feedback 
loops, Proc. National Acad. Sci. USA 104 (2007) 6533-6537. 

[31] R.V. Plykin & A. Zhirov, Some problems of attractors of dynamical systems. Topol- 
ogy Appl. 54 (1993), 19^6 

[32] D. Ruelle, Strange attractors as a mathematical explanation of turbulence. Statisti- 
cal models and turbulence (Proc. Sympos., Univ. California, La JoUa, Calif., 1971). 
Lecture Notes in Physics vol. 12. Berlin: Springer- Verlag (1972). 

[33] D. Ruelle, "Turbulence, strange attractors, and chaos." World Scientific Series on 
Nonlinear Science. Series A: Monographs and Treatises, 16. River Edge, NJ: World 
Scientific PubUshing Co. (1995) 

[34] D. Ruelle, Small random perturbations of dynamical systems and the definition of 
attractors. Comm. Math. Phys. 82 (1981-82), 137-151. 

[35] D. Ruelle & F. Takens, On the nature of turbulence. Comm. Math. Phys. 20 (1971), 
167-192. 

[36] D. Ruelle & F. Takens, Note concerning our paper: "On the nature of turbulence". 
Comm. Math. Phys. 23 (1971), 343-344 

[37] G. Sell, Periodic solutions and asymptotic stability, J. Differential Equations 2 (1966), 
143-157 

[38] G. Sell & Y. You, "Dynamics of Evolutionary Equations." Applied Mathematicatical 
Sciences, 143. Springer- Verlag, New York, 2002. 

[39] S. Smale, Differentiable dynamical systems. Bull. Amer. Math. Soc. 73 (1967), 747- 
817. 

[40] H. L. Smith and H. R. Thieme, Convergence for strongly ordered preserving semi- 
flows, SIAM J. Math. Anal., 22 (1991), 1081-1101. 

[41] E. H. Snoussi, Necessary conditions for multistationarity and stable periodicity, J. 
Biol. Syst. 6 (1998) 3-9. 

[42] E.D. Sontag. Monotone and near-monotone biochemical networks. Systems and Syn- 
thetic Biology 1 (2007) 59-87. 



21 



[43] E.D. Sontag, A. Veliz-Cuba, R. Laubenbacher, and A.S. Jarrah, The effect of negative 
feedback loops on the dynamics of boolean networks, submitted, 2007. 

[44] P.R. Stein & S.M. Ulam, Quadratic transformations. Part I. Los Alamos Report 
LA2305. Washington: Office of Technical Services, U.S. Department of Commerce 
(1959). 

[45] P.R. Stein & S.M. Ulam, Non-linear transformation studies on electronic computers. 
Rozprawy Matematyczne vol. xxxix. Warsaw: Pannstwowe Wydawnictwo Naukowe 
(1964). 

[46] R. Thom, "Stabilite structurelle et morphogenese; essai d'une theorie generale des 
modeles." Reading, Mass.: W. A. Benjamin (1972) 

[47] R. Thom, "Structural stability and morphogenesis; an outline of a general theory of 
models." Reading, Mass.: W. A. Benjamin (1975). 

[48] R. Thomas, On the relation between the logical structure of systems and their ability 
to generate multiple steady states or sustained oscillations. Springer Ser. Synergetics, 
9, (1981) 180-193. 

[49] A.M. Turing, The chemical basis of morphogenesis, Phil. Trans. Royal Soc. B 237 
(1952), 37-72 

[50] T. Ura, Sur les courbes definies par les equations difFerentielles dans I'espace a m 
dimensions, Ann. Sci. Ecole Norm. Sup. (3) 70 (1953), 287-360 

[51] L. Ward, Partially ordered topological spaces, Proc. Amer. Math. Soc. 5 (1954) 144- 
161 

[52] R.F. Williams, Expanding attractors, CoUoque de Topologie Differentielle (Mont- 
Aigoual, 1969), 79-89. MontpeUier: Universite de MontpeUier 



22 



